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We investigate the properties of a degenerate dilute gas 
of neutral fermionic particles in a harmonic trap that inter- 
act via dipole-dipole forces. We employ the semiclassical 
Thomas-Fermi method and discuss the Dirac correction to 
the interaction energy. A nearly analytic as well as an exact 
numerical minimization of the Thomas-Fermi-Dirac energy 
functional are performed in order to obtain the density distri- 
bution. We determine the stability of the system as a function 
of the interaction strength, the particle number, and the trap 
geometry. We find that there are interaction strengths and 
particle numbers for which the gas cannot be trapped stably 
in a spherically symmetric trap, but both prolate and oblate 
traps will work successfully. 

PACS numbers: 03.75.Fi, 05.30.Jp 



I. INTRODUCTION 

The experimental achievement of quantum degener- 
acy in a dilute trapped gas of cold fermionic atoms |l]] 
has stimulated theoretical interest in the properties of 
this fundamental system. Attention has focused on the 
critical temperature || and the detection Q of Cooper 
pairing as well as on the properties of mixtures of vari- 
ous fermionic and bosonic species @|||. Another impor- 
tant problem concerns the interactions between ultracold 
fermions Owing to the exclusion principle, spin- 

polarized Fermi atoms do not interact via s-wave col- 
lisions, whereas they dominate in the low-energy regime 
for bosons and have pronounced effects on the statics and 
dynamical properties of cold boson gases ||. Hence, in 
the absence of low-energy collisions other types of forces 
come into play. 

A good candidate is a dipole-dipole interaction be- 
tween atoms or molecules, not analyzed so far in the con- 
text of cold trapped fermions. Some atoms possess per- 
manent magnetic dipole moments of considerable mag- 
nitude (chromium, for instance, has /J, = 6//g)- It was 
also proposed to induce electric dipoles in atoms |p|,|lO|. 
Huge permanent electric dipole moments occur naturally 
in diatomic polar molecules [[llj . The behavior of atomic 
bosonic dipoles in traps has been investigated in p^| , ^3[ , 
which addressed the question of instabilities in the system 
caused by an attractive component of dipolar interac- 
tions. The conclusion drawn was that a large enough pos- 
itive scattering length (providing repulsive interactions) 
can stabilize a system of bosonic dipoles. For strong 



(e.g. molecular) dipoles, when supposedly the scatter- 
ing length can be neglected, it is the trap geometry that 
plays a crucial role [|0| - the system is stable provided 
the trap assures the domination of repulsive interactions 
in the gas. 

A full quantum mechanical description of the system 
of many interacting fermions is of course very complex. 
But the lessons of semiclassical atomic physics can be ap- 



plied, in particular the Thomas- Fermi approach [|14 15 
and its refinements (see, e.g., @). Its success in de- 
scribing static properties of atoms is well known, and we 
note that, in recent years, these methods were also used 
successfully for studying dynamical processes of atoms 
and molecules in superstrong light pulses . 

Our paper is organized as follows: in Section II, the 
Thomas-Fermi model is revisited with an eye on the 
dipole-dipole forces. The Dirac correction to the inter- 
action energy term is discussed and scaling properties 
derived with the help of the virial theorem. In Section 
III the results of approximate (nearly analytic) and nu- 
merical minimizations of the Thomas-Fermi-Dirac energy 
functional are presented. Unexpectedly, we find that the 
system may withstand larger dipolar forces both for very 
flat and for highly elongated traps. 



II. THOMAS-FERMI MODEL 

A. General considerations 

Let's begin by recalling some basic things, mainly col- 
lected from Refs. [ |l6||i8|]l9| ] , with due attention to the 
changed situation: here fully spin-polarized fermions — 
there fermions with no net spin. The spatial one-particle 
density is denoted by n(r), it is normalized to the total 
particle number N, 



N = I (dr) n(r) 



(1) 



where (dr) = dxdydz denotes the volume element. The 
spatial one-particle density matrix rS^ (?';?") and the 
one-particle Wigner function v(r,p) are related by 



(r ;r ) 



(dp) 



r")/h 



(2) 



1 



The spatial and momental one-particle densities are ob- 
tained by integrating v(r,p) over the other variable, 



n(r) = n^(r;r) 
P(P) 



(dp) 



(2ttH) 3 
(dr) 



(3) 



They are needed for the calculation of the kinetic energy, 

P -2 

Skin = J (dp) f^ P (p) , (4) 

and the external potential energy (of the harmonic trap) , 
Strap = I (dr) \Muj 2 [x 2 + y 2 + (fiz) 2 ] n(r) , (5) 



where M is the mass of the atom species considered, u> is 
the (transverse) trap frequency, and (3 is the aspect ratio 
of the cylindrical trap. The trap is spherically symmetric 
for j3 = 1; for j3 < 1, the equipotential surfaces are prolate 
("cigar shaped") ellipsoids; for p > 1, they are oblate 
( "lentil shaped" ) ellipsoids. 

For the dipole-dipole interaction energy, E dd , we need 
(the diagonal part of) the two-particle density matrix 

(2)1-1 ->' -I I -"\ 

n ( , {r 1 ,r 2 ;r 1 ,r 2 ), 
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where 



- (dr)(dr")V dd (r - r")n (2) (r ,r";r',r") , (6) 
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We note that the contact term, proportional to S(r), is 
required by the condition that the magnetic field made by 
the point dipole be divergence- free |E0|. An alternative 
way of presenting V dd is 



V dd (r) = 



.1 



VV- -4nlS(r) 



(8) 



which is a particularly convenient starting point for eval- 
uating the Fourier transform 



(dr) e ifc -V dd (f) = ^. 



Mo. 



kk 

47Tw — 4-7T1 
K z 



A* ■ (9) 



The vanishing divergence just mentioned is here imme- 
diately recognized, inasmuch as k ■ [• • •] =0. 



B. Thomas- Fermi-Dirac functionals 

The semiclassical approximation now employed — in 
the spirit of what the TFD trio (Thomas @ , Fermi ^ , 
and Dirac pi]]) did, although in a technically different 



manner — is two-fold: We approximate by products 
of factors (Dirac), 

(2) /-»/ — / —II -II \ (l)l-i -a \ (l)l-i -a \ 
n y >(r 1 ,r 2 ]r 1 ,r 2 ) = n y '{r 1 ;r 1 )n y >(r 2 ;r 2 ) 

-nMfcfl)nV><? 3 -fl), (10) 

and 7?/ 1 ) by a brutally simple Wigner function (Thomas 
and Fermi), 



u(r,p) = r](h 2 [6ir 2 n(r)] 3 —p 2 ) 



(11) 



where r)( ) is Heaviside's unit step function. This gives 
the density functional of the kinetic energy as 



E kin [n] = / (dr) 



1 
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[^ 2 n(r)}\ (12) 



and the density functional for the potential energy in the 
trap is E tlap of (|). 

The dipole-dipole interaction energy consists of two 
parts, corresponding to the two summands in (flQ), 



E dd [n]=E^[n]+E^[n] 
with the direct term 



(13) 



4d r) N = \ f (d?)(df)n(r')V dd (?' - r")n(r") (14) 



and the exchange term 



E { d f[n] = -- / (dr')(dr")n (1) (?';r> (1) (?";^') 
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->"i (l) i-i -n\ (l) i— n -<\ 
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xV dd (7'-r") 

= -\ /(dr) j{ds)VM 

xn (1) (r + §s;r- ±s)n (1) (r- ±s;r + |s) . 

(15) 

Now we note that 

n (1) (r + \s;f- ±s)n (1) (r- \s\r + \s) 

'- V {r,p') V {r,f)SP' ~P"^°l h 

(16) 
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depends only on the length s — Is of vector s, not on 
its direction s/s, because — in the TF approximation 
( pd| ) — the product v{r,p')u{j,p") involves only p' = |p'| 
and p" = |p"|. As a consequence, it is permissible to 
replace, in (Jl5|), Vdd (3) by its average over the solid angle 
associated with s, 



— jj-M *(s) 



(17) 



We thus arrive at 
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E dd W = 2 / ( dr ) t 71 ^)] 4^ T M 
= ^ y (dr)(drV(r)^ 
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and accordingly 
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Edd[n] = - / (dr)(dr ) n(r)y dd (r - r )n(r ) 



with 



V dd (r 
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which is Vdd(?) of (j?]) with the contact term removed. 
In (|8|) and (|9|) this corresponds to multiplying the unit 
dyadic T by i @. 

The rotational symmetry of the trap potential in (|5|) 
distinguishes the z axis, and we take for granted that this 
is also the direction of spin polarization, 



Then 



pL = fie z 



T7 MO 2 1 -3( Z / r ) 2 

v **<r) = — — 



(21) 



(22) 



and the whole system is invariant under rotations around 
the z axis. 

The total TFD energy functional is then given by the 
sum of the kinetic energy ( |l2] ) , the potential energy in the 
trap (||), and the dipole-dipole interaction energy (|l9|), 



(dr) 



M 20n' 



■ [6TT 2 n{r)] 3 + \Mw 2 r 2 n{r) 



+- / (dr)(dr)n(r)Vdd(r - f)n(f) . (23) 
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The density that minimizes £'( TFD ) under the constraint 
(El) obeys the nonlinear integral equation 



^[6A(r)]' + iMa, 2 [, 2 +y 2 + (^f] 

+ / (dr)Vdd(r ~ r)n(f') = \Muj 2 R 2 , (24) 



where \Muj 2 R 2 is a convenient way of writing the La- 
grange multiplier for the constraint. 



C. Virial theorems 

Scaling transformations of the form 

n(r) -> A 3+a n(Ar) , N -> X a N 



are consistent with the constraint (|l]). They affect the 
various pieces of E^ TFD ^ in accordance with 



(26) 



n(r ) , 
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- A 3+2a E dd , 




so that (E = E( TFD )) 




(19) 


E = -Ekin + E trap 


+ Edd 




-> A 2 +i Q E kin + 


A 2+a £"trap + 



dd 



(27) 



In the infinitesimal vicinity of A = 1, all first-order 
changes of E originate in the explicit change of N, 
5N = 5\ aN, and therefore 



dE_ 

dN 



aN ^Tr = (2 + H^kin + (-2 + a)E trap + (3 + 2a)E dd 



(28) 



must hold irrespective of the value of parameter a. In 
view of the linear a dependence, we get two independent 
statements, 



a = 



a 



2-E'km — 2E trap + 3E dd — , 

Eki„ + 7£ trap = 3N^ . (29) 



They enable us to express -Ekin, E trap , and E dd in terms 
of E and NdE/dN, 
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E kin = —E - —N 



E trap = --E+-iV 



dN 
3, 8E 



Edd = -8E + 6N 



2 dN 
dE 



dN 



(30) 



In conjunction with (/i = |/i|) 
I x ~q^E — 2E dd , 
w^jE = 2E trap , 



(31) 



they imply that the ground state energy E(M, to, /i, N) 
is of the form (we leave the (3 dependence implicit) 



with 



E(M, to, fx, N) = hujN^e(N^e) 



(32) 



(33) 



Clearly, the dimcnsionlcss number e measures the relative 
(25) strength of the dipole-dipole interaction. The universal 
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function e( ) is to be found numerically (for each /3 value 
of interest). For e = 0, we can solve (f24j) immediately, 

'" (Mio/h) 3 [R 2 ~x 2 -y 2 -(/3z) 2 ] f , (34) 



n(r) 



6tt 2 



where, by convention, [• • •] 2 =0 for negative arguments, 
and 

r = (4&(3Ny (Muj/ny^ (35) 

as a consequence of and so we find 



III. RESULTS 

Owing to its intrinsic semiclassical approximations, 
one expects a few-percent deviation of the TFD energy 
from the true ground-state energy. It is, therefore, not 
really necessary to solve the nonlinear integral equation 
(p!3|). A reasonable variational estimate, in conjunction 
with a few full-blown numerical solutions for comparison, 
will do. 



e = e(0) = -(6/3)3 = 1.363/33 



(36) 



For the spherically symmetric case = 1, E^d vanishes 
in first-order perturbation theory, so that 

/3 = 1: e{Ni e) = 2~h^ +e 2 Nh 2 + 0(Nh 3 ) (37) 
with e2 < 0. 

D. Dimensionless variables 

These scaling laws invite the use of correspondingly 
chosen dimensionless variables, such as 



x = r/a with a = N 1 



h 

Mu 



(38) 



for the position and 

a 3 N 

= TF n ( aS ) or n ( ? ) = s9^/a) (39) 
N a 6 

for the density. The constraint (p]) then appears as 



(dx) g(x) = 1 , 



and the TFD energy acquires the form 

+ \ j (d2) (x 2 + x 2 +(3 2 x 2 )g(x) 

1 ,ri /"/ ,-./■« 3 COS 2 I 

+ -N°e (dx)(dx)g(x) — — ^7— 



(40) 



L J \X — X \" 

(41) 

where 9 is the angle between the polarization direction 
(the #3 direction) and the relative position vector x — x' , 

(x — x) ■ Jl = \x — x'\[xcos9 = (x% — x' 3 )n . (42) 

The minimum of the scaled TFD energy is e(Nee) of 
( |32| ) (with its implicit /3 dependence); it is obtained for 
the g(x) that obeys the dimensionless analog of (|24|), 



\[&n 2 g{x)Y + \(x\+x 2 2 +(3 2 x 2 ) 

» T i f, 1 — 3 cos 2 9 , „ . 

+me / (dx') ^ g(x') = \X 2 , (43) 



where the value of X — R/a is determined by the con- 
straint (Eoh. 



A. Gaussian variational ansatz 



The Gaussian ansatz 



(2ir) 2 e 2 



\k 2 (x\ +xl +7 2 x§) 



(44) 



is convenient. In addition to its aspect ratio 7, the shape 
parameter, it contains the scale parameter k, so that the 
virial theorems of Sec. II C will be obeyed for the optimal 
choice of k. For this scaled density, the scaled kinetic 
energy is given by 



and the scaled value of -Etrap 



= 2 3 3 e 5 2 tt 3 K z j 3 
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(45) 



(46) 



They exhibit the anticipated dependence on the scale pa- 
rameter k, and so does the dipole-dipole interaction en- 
ergy, 

-1 



E dd 
hujNi 



4^¥ 



c 2 -c 4 



1 - C 2 + 7 2 C 2 



(47) 



This integral can, of course, be evaluated in terms of 
elementary functions, but as it stands we see immediately 
that the integrand, and thus the integral, is positive for 
< 7 < 00, so that 

-Edd > for 7 > 1 (oblate Gaussian), 

Edd = for 7=1 (spherical Gaussian), 

i?dd < for 7 < 1 (prolate Gaussian), (48) 

consistent with the expectation that an oblate density 
has a larger magnetic interaction energy than a prolate 
density. More explicitly, then, we have 



E, 



dd 



huNi 4^ 



Ne£J 
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sin 2 1? 3 ' 
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> 1, 



< 1 



(49) 
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FIG. 1. Stability diagram. The dots, connected by a solid 
line to guide the eye, define the border between system param- 
eters (/3: aspect ratio of the trap; Nee: effective interaction 
strength) for which i?( TFD > is bounded from below (stable) 
or not bounded (unstable). The inset shows a blow-up of the 
region of pronouncedly prolate traps (j3 < 0.5). For oblate 
traps with (3 = 5.2 or larger, the system is always stable, ir- 
respective of the value of Nee. This figure presents results 
obtained in the Gaussian approximation, and so do all other 
figures. 



The Gaussian density ( fl4|) cannot mimic the e = so- 
lution ( p4| ) very well. But nevertheless, the resulting es- 
timate of the e — energy, 

Gaussian: e = 2~53if 5-1^5/35 = 1.42/3* (50) 

is only 4.4% in excess of the correct value (|36|). This 
accuracy is sufficient for our purposes; and, in any case, 
an error of a few percent is a small price for the enormous 
simplification that the Gaussian ansatz brings about. 

Dipolar interactions are partially attractive and par- 
tially repulsive, depending on the configuration of the 
dipoles. One should keep in mind the simple situation of 
two dipoles in the plane perpendicular to their polariza- 
tion, which repel each other, as opposed to the situation 
of two attracting dipolar particles placed along the di- 
rection of their polarization. Extending this picture to a 
cloud of trapped dipoles one would expect that attrac- 
tion dominates in prolate traps, and repulsion in oblate 
traps (provided the dipoles are polarized along the trap 
axis as is the case here). In the case of predominant 
attraction one may surmise that instabilities occur. By 
varying the two system parameters (Nee and f3) we have 
investigated the issue of stability, see Fig. 0. 

From this stability diagram one concludes that, for 
oblate traps ((3 > 1), the bigger the trap aspect ratio, 
the bigger values of the dipole parameter Nee can be 
stabilized. In fact, we found numerically that fermions 
form stable configurations for (3 > 5.2 irrespective of the 
strength of their dipole interaction (an analogous effect 
was observed for dipolar bosons by Santos et al. [pi). 
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FIG. 2. Dependence of cloud parameters n (cloud size; top) 
and 7 (cloud shape; bottom) on f3, the aspect ratio of the trap. 
The lines refer to different values of the interaction strength: 
Nee = (solid lines), Nee = 1 (dashed lines), Nee = 2 
(dash-dotted lines), Nee = 2.4 (dotted lines). 



On the other hand, one might naively expect an exactly 
opposite effect in prolate traps (/? < 1) where attractive 
interactions should dominate. This is not quite true - 
indeed for moderate trap aspect ratios (1 > (3 > 0.5) the 
critical value of the dipole parameter is smaller, but in 
traps that are very soft (/? < 0.5) in the z direction of ro- 
tational symmetry we observe an increase of the critical 
value of N^e (see the inset in Fig. [j]). This can be under- 
stood with the help of the following argument. We note 
that the dipole-dipole energy term vanishes for a uniform 
density distribution. As the trap is made softer in the po- 
larization direction, the shape of the cloud along the soft 
axis becomes more and more uniform, contributing to 
the interaction energy to a lesser extent (this argument 
was also used to interpret our earlier results for bosons 
interacting via contact and dipole-dipole forces, see |l3| ). 

The dependence of 7 and k on (3, shown in Fig. Mis 
consistent with this argument. We see that 7 decreases 
with decreasing whereas k increases. Accordingly, the 
cloud gets stretched along the z axis of symmetry, and the 
diameter of the circular cross section in the x, y plane (oc 
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FIG. 3. Aspect ratio 7 of the cloud as a function of the in- 
teraction strength N'e'e, for various values of the aspect ratio 
P of the trap. Note that 7 = /3 for TV 5 e = 0. The dashed line 
is for P = 5.2, for which the cloud is stable for all values of 



k^" 1 ) is reduced. At the center of the cloud, we thus have 
a relatively large volume of (almost) constant density, 
and the inhomogeneous parts of the cloud are relatively 
far apart. Taken together, these geometric features lead 
to a rather small dipole-dipole interaction energy. 

Whereas the stability of the system considered can be 
well understood, the spatial behavior of the fermionic 
cloud, especially near the collapse, seems to be much less 
intuitive. As we approach the critical parameter values, 
the aspect ratio 7 of the cloud decreases and the cloud 
becomes elongated in the attractive z direction. This 
type of behavior is general in the sense that it does not 
depend on the trap aspect ratio, see Fig. ||. It is only for 
traps with /3 > 5.2 that 7 reaches an asymptotic value 
(dependent on 0) for extremely large values of Nee. 

The dependence of the dipolar energy Edd on the 
dipole parameter Nee and the trap geometry is also of in- 
terest as it is the quantity responsible for the (in)stability 
of the system, see Fig. 0. For all prolate traps, Edd re- 
mains negative approaching some critical value at the 
collapse point. For (3 < 5.2, the dipolar energy can be 
positive for moderate dipole parameters, but if their val- 
ues are large enough, Edd turns negative indicating the 
collapse. For (3 > 5.2, Edd is always positive and in- 
creases as a function of Nee. 

Finally, we take a look at the total TFD energy, see 
Fig. ||. For (3 < 5.2, the system becomes unstable at the 
critical value of Nee. Consistent with Fig. El we observe 
that larger values of Nee are supported for f3 — 0.1 and 
f3 = 1 than for (3 = 0.5. 

Let us now discuss the dipole parameter N^e and give 
some typical values for it. Owing to the A~ dependence, 
one can locate the experimental system in various regions 
of the stability diagram of Fig. [I] not only by choosing (or 
inducing, as proposed for bosons f|,[l(|) a specific value of 



/i, but also, to some extent, by varying N. One could also 
exploit the uj dependence of e cx y/w, which is particularly 
relevant for optical traps with tight confinement ]22] j. 

For the parameters of the fermionic chromium isotope, 
lu = 300 Hz and N = 10 6 , one obtains N*e = 0.012, 
which is a very small value. Therefore, the conclusion 
for atoms possessing even relatively large magnetic dipole 
moments is that irrespective of their number and the 
trap frequency they will always remain stable against 
collapse. The following amusing analogy offers a good 
reason for this observation [|| . Let us compare the char- 
acteristic sizes of the noninteracting Fermi gas and the 
Bose-condensed atomic gas interacting via a repulsive 
contact potential. In the Thomas-Fermi approximation 
the appropriate quantities, in units of yjh/ (Mu>), read: 

^Formi = (48A r )s for fermions and -Rbosc = (l5Na )i for 
bosons , where a is the scattering length. By equating 
the two sizes one can calculate the effective, A~ depen- 
dent, scattering length due to the exclusion- principle- 
induced repulsion: ao ~ 1.68A~e. For typical num- 
bers of atoms, N = 10 3 • • • 10 6 , this ao is huge on the 
scale set by typical scattering lengths for bosonic atoms 
(ao ~ 10~ 3 ). Once we realize how strong is the repul- 
sion that originates in the Fermi statistics, we under- 
stand that small atomic magnetic dipoles can hardly have 
a noticeable effect on the behavior of dipolar fermionic 
gases. However, for polar molecules the situation is dif- 
ferent. For a trap frequency of uj = 300 Hz and N = 10 6 
molecules of mass m w 100 a.m.u. (the typical mass of 
an alkaline dimer) and a typical electric dipole moment 
of 1 Debye, one reaches Nee ss 11.5, which may very 
well put the system into the unstable regime - see Fig. [l]. 
In this situation, a sufficiently large (3 value will stabilize 
the system. 

In order to assess the quality of our variational results 
we have computed exact numerical solutions of Eq. . 
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FIG. 4. Dipole-dipole interaction energy E^a as a function 
of the interaction strength Nee. The solid lines are for trap 
aspect ratios (3 — 0.5, 1, 1.5, . . . , 4.5; the dashed line is for 
fi = 5.2. 
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FIG. 5. Normalized TFD energy as a function of the inter- 
action strength for various trap shapes. The universal func- 
tion e(N^e) of (p^), normalized to its initial value eo = e(0), 
is shown for trap aspect ratios f3 — 0.1, 0.5, 1, 1.5, . . . , 4 (solid 
lines) and f3 — 5.2 (dashed line). Note the horizontal slope of 
the /3 = 1 line (spherical trap), as required by (p7|). 



This equation was solved for the density distribution g(x) 
sclf-consistently starting from the known analytical re- 
(|34|) for a non-interacting (e = 0) Fermi gas in a 
and slowly increasing the dipole parameter Nee. 
For each value of N^e the solution was iterated until 
convergence was reached. Then, the value of the dipole 
parameter was slightly increased. In order to compute 
the dipole (integral) term we note that it has the form 
of a convolution. Thus, it can be conveniently evaluated 
in the Fourier space where it is a simple product of the 
Fourier transforms of the density (computed numerically 
with the aid of an FFT) and the interaction potential, 
the latter being known analytically ]r^] : 




(dx 



^ p 1 ? ' x 



1 - 3cos 2 6> 



4?r 



(1-3 cos 2 a), (51) 



where a is the angle between the Fourier variable q and 
the z direction. In order to assure that the integral term 
is evaluated accurately we used a Gaussian distribution 
for comparison and chose the grid parameters accord- 
ingly. 

Our numerical calculations, performed in three dimen- 
sions, were quite demanding so we limited their use to 
a check of the main features of the stability diagram in 
Fig . [|. The solutions obtained satisfy the virial relations 
(E9) very well, and total energies obtained numerically 
are always below the corresponding values from the vari- 
ational analysis. The critical value of the dipole param- 
eter for a spherical trap is Nee = 1.96 as compared to 
2.55 obtained by the variational calculation. We also 
confirmed the effect of increase of the critical interaction 
strength that we found, in the Gaussian approximation, 
for prolate traps: for (3 = 0.07 the critical value of Nee 



is 2.75(> 1.96) as compared to 2.81(> 2.55) obtained an- 
alytically in the Gaussian approximation. 

Now two remarks about possible extensions of our work 
are in order. Firstly, the results presented in this pa- 
per describe the situation of a dipolar fermionic gas at 
the temperature T = 0. An interesting subject of study 
would be extension of our theory to finite temperatures. 
Secondly, there exists a parallel approach in the Thomas- 
Fermi model, namely the one in the momentum space 
p8|,^9[ . As many experiments with cold gases yield their 
momental characteristics, investigation of this alternative 
approach also presents an attractive theoretical task. 
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